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In weakly nonlinear dispersive systems, 
solitons are spatially localized solutions which 
propagate without changing shape through a deli¬ 
cate balance between dispersion and self-focusing 
nonlinear effect^. These states have been ex¬ 
tensively studied in Bose-Einstein condensates 
(BECs), where interatomic interactions give rise 
to such nonlinearitie^. Previous experimental 
work with matter wave solitons has been limited 
to static intensity profiles. The creation of matter 
wave breathers —dispersionless soliton-like states 
with collective oscillation frequencies driven by 
attractive mean-field interactions—have been of 
theoretical interest!^ due to the exotic behaviour 
of interacting matter wave systems. Here, using 
an attractively interacting Bose-Einstein conden¬ 
sate, we present the first observation of matter 
wave breathers. A comparison between experi¬ 
mental data and a cubic-quintic Gross-Pitaevskii 
equation (GPE) suggests that previously un¬ 
observed three-body interactions may play an 
important role in this system. The observation 
of long lived stable breathers in an attractively 
interacting matter wave system indicates that 
there is a wide range of previously unobserved, 
but theoretically predicted, effects that are now 
experimentally accessibld^. A compelling exam¬ 
ple is the recently predicted quantised tunnelling 
of a breather through a barrieiP. 

The study of solitons has found application in fields 
as diverse as nonlinear optic^, magnetic material^, and 
atom interferometr}El. Both darli^ and bright matter 
wave solitons have been realized experimentally and have 
been shown in exotic contexts such as colliding solitonJ^, 
soliton formation during cloud collaps^I^, and the colli¬ 
sion of bright solitons with a barriei^i^. 

Breathers are thought to be a similarly ubiquitous 
phenomenorP31. They have been experimentally realised 
in a variety of systems including Josephson-junctio n lat- 
tices, spin lattices, and optical waveguide array^^^^. 
Breathers are also thought to play a role in more exotic 
contexts such as DNA denaturatiorP^. While breathers 
have been predicted to exist in matter wave system^, 
they have not been observed until now. 

To create the matter wave breathers, an attractively 
interacting Bose-condensed cloud of ^^Rb atoms is made 


to undergo collective quadrupole oscillations inside an 
anisotropic cigar-shaped trap (see Figure [^ . A se¬ 
quence of measurements are made on the excitations of 
trapped clouds at varying s-wave scattering lengths, a^, 
the critical parameter determining the lowest order in¬ 
teractions of the atoms. In the regime where < 0, sta¬ 
ble quadrupole oscillations above the frequency of those 
measured at = 0 are observed, indicating that the 
quadrupole mode in this regime is being driven by the 
attractive s-wave interactions, confirming the existence 
of breather 

In addition to the observation of matter wave 
breathers, we find quantitative evidence of physics not 
captured by two-body scattering alone. Firstly, a mea¬ 
surement of the quadrupole oscillation frequency against 
the fundamental centre of mass (CoM) mode of the trap 
in the absence of two-body interactions {as = 0) gives an 
unexpected, but previously predicted ratio. Secondly, at 
negative s-wave scattering lengths (a^ < 0) the trapped 
condensate is also observed to be stable over a larger 
range of scattering lengths than previously predicted. 



FIG. 1. a) The UHV apparatus, with the glass cell centred. 
Solenoids top and bottom provide a magnetic bias field of 
~155G while intersecting laser beams at 1064nm generate a 
potential that traps the atoms, b) Schematic of a BEG un¬ 
dergoing a quadrupole oscillation inside a cigar-shaped op¬ 
tical dipole trap. The quadrupole oscillation occurs around 
the weaker axial frequency of the harmonic optical trapping 
potential holding the atoms, c) The quadrupole collective 
normal mode. Red and blue isosurfaces are out of phase con¬ 
tributions to the condensate density oscillations. The mode 
was tomographically reconstructed from dual axis imaging 
data using principal component analysis and an inverse radon 
transform with imposed cylindrical symmetry. 
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Finally, an anomalous phase shift is observed between 
the axial and radial quadrupole oscillations in expan¬ 
sion data. It is shown that the presence of a repulsive 
three-body interaction quantitatively explains the obser¬ 
vations. 


Initially, the quadrupole mode oscillation was probed 
in the absence of two-body interactions (a^ = 0). The 
scheme described in Methods was used to excite the 
quadrupole mode (ooq) at = 0 and the oscillation 
frequency was measured by fitting the averaged widths 
to a decaying sine curve. This was compared to a di¬ 
rect measurement of the axial trapping frequency (oOz) 
made by using Bragg transitions to excite the centre of 
mass (CoM) mode. Data was acquired by alternating be¬ 
tween quadrupole and CoM measurements to reduce the 
chance of systematic drift influencing the measured ratio 
It was found that at = 0, ujq = (1.67±0.03)ci;;2 
with ujz = 27r x (7.55 ± 0.15)Hz. This ratio differs sig¬ 
nificantly from the noninteracting gas {as = 0) result 
which predicts a ratio of 2 between quadrupole and CoM 
modes. Measurements made at positive scattering length 
for both the ^^Rb BEC and a ^^Rb BEC in the same 
trap give results consistent with theoretical predictions 
{ujQ = ^/K/2uJzj^• ^^Rb is observed to have a significant 
three-body loss around the Eeshbach resonance at 155GP^ 
and, intriguingly, it is predicted that there should also be 
significant three-body scatterin^^. The ratio measured 
here for the noninteracting gas is in excellent agreement 
with a theoretically predicted frequency ratio which ac¬ 
counts for non-zero three-body interactions at 
(see supplementary material). 


Having investigated the = 0 case, the quadrupole 
mode oscillation was excited at varying negative scatter¬ 
ing lengths, as < 0. In these experiments, the trapped 
cloud of ^10^ atoms is observed to be stable against 
condensate collapse over a broader range than expected 
from a mean field theory with only two-body scattering 
(see supplementary material for discussion and data). In 
this regime, the negative s-wave scattering length will 
significantly modify the collective excitation spectrum. 
Eigure shows the axial and radial oscillations of an 
N=1.5 X 10^ BEC, induced by jumping the scattering 
length from = 40ao to = — 13ao, where ao is the 
Bohr radius. The oscillations are stable for over 400 ms. 
Eor as = — 13ao the quadrupole mode frequency is de¬ 
termined to be ujQ = 27r X (11.54 ±0.12) Hz. A small 
phase shift of (0.24 ± 0.06) tt between the axial and ra¬ 
dial widths is present. This shift arises during the 20ms 
of free space expansion before absorption imaging of the 
cloud. In the supplementary material, it is shown that 
this shift during the expansion is quantitatively predicted 
by an extended GPE containing a quintic ternP^ 



FIG. 2. The measured radial (r) and axial (z) widths of the 
oscillating cloud in the cross trap after the scattering length 
was jumped to — 13ao. Each point is a weighted average of 
multiple shots (~4). Fitted decaying sine curves with ujq = 
27r X (11.54 ± 0.12) Hz are also plotted. Inset are absorption 
images at 77 ms, 92 ms and 122 ms showing the cloud changing 
from maximum to minimal axial extension. The radial and 
axial widths oscillate almost out of phase, with a small phase 
shift of (0.24 ± 0.06) tt due the free space expansion of the 
condensate prior to imaging (see supplementary material). 

= {-^^^ + V{r)+g2N\^{r,t)\^ 

where r = (x, z) is a position vector, N is the con¬ 
densate population, m is the mass of an ^^Rb atom, V{r) 
is the trapping potential, g 2 = is the two-body in¬ 

teraction term, and gs is the three-body interaction term. 
For a cylindrical cigar-shaped trap, as in our system, the 
potential is given by V{r) = ^ ± 0Jr{x^ ± y^))- 

By conducting measurements from = 0 through a 
range of negative scattering lengths the dependence of 
the quadrupole frequency on can be determined (Fig¬ 
ure]^. The results from two separate experiments (the 
optical trap was rebuilt between data sets) using con¬ 
densed populations of = 1.5 x 10^ and = 1 x 10^ 
are plotted. The data for different were taken out of 
order to remove systematic effects on the measured fre¬ 
quencies. As as becomes more negative the oscillation 
frequency becomes greater than the measured interac¬ 
tionless {as = 0) quadrupole mode. An observation of a 
collective oscillation frequency higher than the measured 
interactionless case is the hallmark of a breather. 

An analytic expression for the quadrupole excitation 
frequency for a cubic-quintic GPE has previously been 
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• Set 1 ■ Set 2 m Analytic □ Simulation 



FIG. 3. The measured oscillation frequency against the scat¬ 
tering length for set 1 (iV = 1.5 x 10^) and set 2 (iV = 1 x 10^). 
Frequencies are normalised by ujz which is calculated from 
the interactionless (as=0) oscillations. Dashed lines indicate 
the theoretical frequencies of equation § using parameters 
{N=1.5 X 10^, ujz = 27 t X 5.88 Hz, cjr = 27r x 77 Hz, Re[^ 3 ] = 
1.23 X and {N=1.0 x 10"^, = 27r x 6.61 Hz, 

tJr = 27r X 77 Hz, Re[g 3 ] = 3.3 x 10”^® ?icm®s“^} for sets 1 
and 2 respectively. Open data points are simulated frequen¬ 
cies for the same parameters. The inset figure demonstrates 
the oscillation measurement taken for each data point. 


derived^ 


luq — ^2 + ^kjQ^2 + — 2P ^ 2,3 

~‘^\J ^ 2/c J6,2 — + P'^ 2 , 3 ) + 8 (PJ3^2 + 2/c J5^3) 

( 2 ) 


where A = is the trap aspect ratio, P = 

y^2/TTNas dimensionless two-body interaction pa¬ 


rameter, and k = ^6 i^ ^ dimensionless three- 

body interaction parameter. Jij = 1 / {AuI^qU^q) , where, 
respectively, Uro and Uzo are the equilibrium axial and 
radial widths of the cloud in dimensionless units. 

As highlighted in the supplementary material, a num¬ 
ber of independent quantitative pieces of evidence point 
to a three-body interaction parameter playing a role in 
the system dynamics. In the following, the experimen¬ 
tal data from Figure are analysed assuming that there 
is indeed such a three-body interaction. To compare 
with theory, equation <§ was fitted to the experimen¬ 
tal data by varying the value of ^ 3 , with ujz being set by 
the measured interactionless oscillation. Figure shows 
the resulting fit for the two sets of data. Fitted values 
are Re[ 5 f 3 ] = (1.23 ± 0.09 ± 0.2) x 10“^^ hcm^s~^ and 
(3.3 ± 0.4 ± 1.0) X 10“^^ hcm.^s~^ for sets 1 and 2 respec¬ 
tively, well within the range predicted by theor}i^^M^ (see 
further discussion in the supplementary materials). The 


first uncertainty quoted for Re[^ 3 ] is statistical while the 
second is from systematic effects such as atom number 
variation. Quadrupole oscillation frequencies obtained 
from numerical simulations of Equation are also pre¬ 
sented in Figure showing excellent agreement with the 
analytical model. For the data in Figure the ratio 
predicted by equation at = 0 is ujq/uJz = 1-8 for 
both sets. Equation!^ also predicts the measured ratio of 
CoM to quadrupole oscillations from measurements with 
Bragg {(jJqIuJz = 1-7 for the experimental conditions in 
that measurement). 

All axial oscillations in the range = [—30ao, Oao] are 
seen to occur around a stable mean width (up to a small 
decay from atom loss), as in Eigure[^ It is also noted 
that as Us decreases the quadrupole mode frequency be¬ 
comes higher than the = 0 case indicating that the 
oscillations are driven by nonlinear self-interac tion o f the 
cloud rather than the trapping potential alon d^^ * ^^ . 

In conclusion, the first observation of a matter wave 
breather driven by attractive two-body interactions was 
made. The dynamics of these breathers were found to 
agree with the predictions of a cubic-quintic GPE equa¬ 
tion with repulsive three-body interactions. The ex¬ 
tended GPE was found to closely model the stability of 
the condensate at negative scattering length, the mea¬ 
sured values of the breather oscillation frequencies, as 
well as the phase shift due to cloud expansion, yielding 
an estimate for the three-body scattering parameter 
in ^^Rb to be of the order 

It is clear from the data presented here that further 
experiments aimed at elucidating the three-body interac¬ 
tion are required. Although this higher order interaction 
quantitatively models the experiment very well, a key re¬ 
sult indicates that there is more to the story: when the 
two sets of data presented in Eigurej^are used to con¬ 
strain both the analytic and numerical models, different 
values for the g^ parameter are found. The results from 
each set are different by a factor of roughly two and is out¬ 
side our statistical and known systematic uncertainties. 
A re-examination of previous experiments in the attrac¬ 
tive two-body regime for ^^Rb, such as the BosenovgP^ 
and soliton formation?^, may provide additional insights. 
Eurthermore, it may be the case that higher order inter¬ 
actions are observable in other BECs of other species. 
Indeed, evidence of such interactions has recently been 
found in exciton-polariton condensate^. 

Work is ongoing in a number of directions. Eirst, to 
extend this study to freely propagating breather states, 
by removing the axial trapping potential and releasing 
the atoms into an optical waveguid^. The Bragg system 
used to excite the sloshing mode of the cloud can also be 
applied to transfer quantised centre of mass momentum 
to the breather state, as well as splitting the cloud into 
multiple momentum states with varying ratioP^. Com¬ 
bined with potential barriers created by blue-detuned, 
tightly focused lasers, we are now in a position to be- 
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gin studying the exotic nonlinear physics of multiple in¬ 
teracting matter wave solitons. Immediate possibilities 
include examining collisions of excited and ground state 
solitons, the study of breather solitons interacting with 
potential barriers, and the examination of the stability 
and breakup of these new states. 


METHODS 

Experimental 

The experimental apparatus is described in depth in 
previous worlP^I, in summary, a combined 2D and 3D 
MOT system collects and cools both ^^Rb and ^^Rb 
atoms. The atoms are loaded into a magnetic trap and 
undergo RF evaporation before being loaded into an op¬ 
tical cross trap. The cross trap consists of intersecting 
1090 nm and 1064 nm lasers with approximate waists of 
300jam and 250yam (half width at 1/e^ intensity), respec¬ 
tively. After loading, the magnetic trap coils are switched 
from anti-Helmholtz to Helmholtz configuration. This al¬ 
lows the 5-wave scattering length of the cloud, a^, to be 
tuned using a Feshbach resonance. Setting the ^^Rb scat¬ 
tering length to zero (a^ = 0) while ramping down the 
cross trap intensity allows sympathetic cooling to remove 
the remaining ^^Rb atoms while minimising three-body 
recombination losses in ^^Rb. A further period of evapo¬ 
ration with as = 300ao creates a pure ^^Rb BEC with 
atom number of the order of 10^. This is the initial 
condition for all experiments. Two orthogonal absorp¬ 
tion beams, one of which is co-axial with the vertical 3D 
MOT beam, allow the cloud to be imaged after 20ms of 
ballistic expansion (a^ = 0) in free space. 

The method of inducing oscillations in the condensate 
is as follows. The initial BEC is held in the cross trap 
while the scattering length is adiabatically ramped to a 
value above that at which the collective oscillations are 
to be measured (typically ^ 40ao). The scattering length 
is then rapidly jumped (in < lOOps) to the desired value. 
Because the spatial width of the ground state at = 40 
is much larger than that at < 0, the cloud begins to os¬ 
cillate. The trap is then turned off and the cloud allowed 
to expand at = 0 before an absorption image is taken. 
By holding the condensate in the cross-trap for different 
times, the temporal dependence of the cloud width can 
be measured. Multiple runs of the experiment are taken 
for each time step to build up statistics. Eor each absorp¬ 
tion image, a 2D Gaussian is fitted to the atomic profile. 
The radial and axial widths, peak intensity and position 
are extracted from multiple shots at each time step and 
an average is constructed by using the error in each fit 
as weights. 


Simulation 

Simulations were conducted to complement the ana¬ 
lytic approach. Equation 0 is numerically solved using 
a 4th-order Runge-Kutta symmetrised split-step Eourier 
method. The dimension of the simulations is reduced by 
exploiting the cylindrical symmetry of the condensate to 
solve the 3D system in 2D cylindrical co-ordinates by a 
Hankel transfornP^. 
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SUPPLEMENTARY MATERIAL 
I. Extracting oscillation frequencies 

The observation of trap oscillations frequencies for 
tts < 0 in required extreme stability of the appara¬ 
tus. Condensates with atom number variation at the 5% 
level over 12 hours and low noise, multi-directional imag¬ 
ing were required to obtain a strong signal. All data sets 
were collected non-sequentially and different data sets 
were taken on multiple days and under different experi¬ 
mental conditions. The Feshbach magnetic field is char¬ 
acterised to significantly below an uncertainty of 1 mG 
using radio frequency spectroscopy of the trapped cloud, 
resulting in a precise, absolute knowledge of at the sub 
O.lao leveP. The conclusions of the paper were drawn 
from a number of different analysis techniques. Critical 
to all analysis was the simultaneous absorption imaging 
of the atom cloud in two orthogonal directions, one with 
an angle of 11.7° with the axial direction and one per¬ 
pendicular. The orthogonal images obtained allowed the 
identification and measurement of all the lower order os¬ 
cillations of the trapped cloud. The top panels of Figure 
l^sup) shows oscillation data extracted by simply fitting 
a Gaussian profile to the absorption image of the cloud 
after 20ms of expansion at zero scattering length. The 
lower panels of Figure [^sup) shows the outcome of prin¬ 
cipal component analysiJ^PCA) on the absorption im¬ 
ages combined from the two orthogonal directions: one 
vertical imaging both the axial and horizontal radial di¬ 
rections (upper image, axial direction is vertical), the 
other imaging down the waveguide giving a radial cross 
section (lower image). PGA supplements a more tradi¬ 
tional analysis and can recover good approximations to 
the shape and frequency normal modes of a BEG from a 
time-series of absorption images. 

Fitting directly to the cloud and PGA give consistent 
results regarding oscillation frequencies. 


II. Oscilations in the presence of three-body 
interactions 

A variational method has previously been used to cal¬ 
culate the collective excitation spectrum for a BEG in 
an anisotropic trapP^. This method agrees well with ex¬ 
perimental results at positive scattering length^. Eor 
Us < 0 this model predicts rapid increase of the oscil¬ 
lation frequency and then collapse of the cloud before 
as ~ — 2ao (for the trap frequencies and atom numbers 
used in these experiments). However, the condensates 
observed here are stable at negative scattering lengths far 
beyond this prediction, indicating that there is physics 
not captured by this model which contributes to the 
condensates’ stability. Attractive mean-field interactions 
will cause the formation of high density regions during 
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FIG. 4. (sup) This data was obtained as follows: from 
BEG creation a 200ms ramp to 40ao in the optical cross trap. 
The scattering length is then jumped to Oao. The cloud was 
moved away from the center of the trap using two 2Tik Bragg 
pulses 0.7ms apart. Images were then taken in 1ms time steps 
and the cloud shape observed. Top panels show experimental 
data for both the CoM oscillation and the quadrupole mode 
at tts = 0. The ratio of these frequencies is given in the text. 
The lower absorption images show the outcome of PC A. The 
first three principal components (PC) are shown where z is 
axial, y is vertical, and x is horizontal, a) Axial CoM mode: 
vertical oscillations in the upper picture indicate the axial 
CoM mode. The small signal in the lower image is due to the 
small angle the second camera makes with the waveguide, b) 
Monopole Mode: in phase axial and radial oscillations, and c) 
Quadrupole Mode: out of phase axial and radial oscillations. 


collapse where higher order interactions will become im¬ 
portant. At large negative scattering lengths, where an 
attractive two-body interaction would otherwise cause 
the condensate to catastrophically collapse, a repulsive 
three-body interaction will stabilise it. 

The three-body term is a complex number with Re[ 5 f 3 ] 
describing the three body scattering parameter and 
Im[^ 3 ] describing the three-body recombination loss. 
The study of trap loss rates as well as the dynamics of 
cloud collapse have so far focused on Im[^ 3 ] with ex- 
perimen tal va lues ranging from 10“^^ to 10“^^ hcm^s~^ 
for 85j^|j20l34H36] Experimental quantification of Re[^ 3 ], 
however, has largely been neglected. Theoretical 
predict ion^^^M^ have ranged from 10“^^ hcm^s~^ to 

By comparing atom number decay curves from experi¬ 
ments and simulations an estimate of the imaginary part 
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of gs was obtained, Im[ 5 f 3 ] = (6 ± 2 ) x 10 h cm^s 
which is well within the range of published valueJ^^^^^M^. 


III. Oscillations of a noninteracting cloud 

The expected frequency of the quadrupole mode when 
as = 0 and ^3 = 0 is ujq = 2uJz However for 

non-vanishing ^ 3 , both Equation 2 and numerical re¬ 
sults from Equation 1 show that this frequency is low¬ 
ered by as much as 15%, depending on the system pa¬ 
rameters. While the quadrupole mode frequency is re¬ 
duced, simulations showed that centre of mass (CoM) 
oscillations (expected to be at ujz) remained unchanged. 
A repulsive 3-body interaction term at = 0 predicts 
a ratio ujqIuJz < 2. The measured ratio at = 0 is 
ujqIuJz = 1-67 ± 0.03 with uoz = x (7.55 ± 0.15)Hz, is 
consistent with Equation 2, for our system parameters 
and Re[^ 3 ] = 3.3 x 10“^^ In Eigure l(sup) 

the experimental data for the quadrupole mode is shown 
compared to the measurement of the CoM. It has previ¬ 
ously been suggested that the two body interactionless 
regime (a^ = 0 ) is the ideal place to search for additional 
3-body physicThus, this measurement alone con¬ 
stitutes an observation that there is significant 3-body 
scattering in ^^Rb near the Eeshbach resonance at 155G. 
Such a shift has not been obser ved fo r oscillations of a 
Li condensate in the same regimd^^H^. While Li and Rb 
systems are similar, the details of the scattering poten¬ 
tials and atomic structure are critical in determining the 
interactions for any particular atonP^. As such an ab¬ 
sence of three-body scattering in Li does not preclude its 
presence in Rb. It should also be noted that no attempt 
has previously been made to measure the excitation spec¬ 
trum of a non-interacting gas of ^^Rb. 


IV. Stability against collapse 

Perhaps one of the most intriguing indications that 
three-body interactions are playing a role in this exper¬ 
iment is the observation of large stable atom clouds at 
Us < 0 . Stable attractively interacting condensates are 
observed to exist beyond the limit predicted by a mean 
field theory considering only two-body interaction^^. 
Eor parameters of data set 2 in the main text Eigure 
3, collapse is observed at 30ao and indicated by a sig¬ 
nificant atom loss as well as an increase in the thermal 
fraction of atoms in the trap. This collapse point is more 
than a factor of 10 larger than predicted by two-body 
scattering alone. Numerical simulations of Equation I 
indicate that a non-zero Re[^ 3 ] stabilises the atom cloud, 
allowing condensates to exist in the range observed here. 
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FIG. 5. (sup) Upper frame: The maximum predicted stable 
atom number, and number of solitons that would be produced 
in the trap used in the experiment, when Re[^ 3]=0 for neg¬ 
ative scattering lengths. Inset is the expected quadrupole 
oscillation frequency when Re[^ 3]=0 for the system parame¬ 
ters of both sets in Figure 3(main). Lower frame: a simple 
model for coupled ID oscillators indicates that for increasing 
number of particles the breathing frequency decreases. 


The observation of condensates violating the predic¬ 
tion of models incorporating two-body scattering alone 
is not unprecedented. A paper by the JILA group re¬ 
ported condensates with up to 10 times the atom num¬ 
ber predicted by GP theor}®. The interpretation of that 
experiment was that solitons were forming in-trap post 
collapse. The absorption images presented in the pa¬ 
per appeared to support this. Numerical simulation of 
the conditions of that experiment, including our inferred 
Re[ 5 f 3 ] interaction, indicate that violent dynamics occur 
during the 2 ms ballistic expansion, because the cloud is 
directly released from the trap at —400ao. The dynamics 
give rise to structures that appear similar to solitons. 

The top frame of Figure [^sup) illustrates the results 
of applying the JILA ‘multi-solition’ model to the condi¬ 
tions here. At the largest negative scattering length used 
for the breather, II solitons would be required to add up 
to the observed stable atom number. Furthermore, a sim¬ 
ple model of coupled oscillators indicates that a soliton 
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chain breathing mode frequency should dramatically de¬ 
crease with an increasing number of discrete solitons (see 
the lower panel of Figure |^sup)). The observation in this 
paper is that the breather frequency increases smoothly. 

It could be argued that the imaging system used here 
is simply not able to resolve soliton trains. In order to in¬ 
vestigate this hypothesis, further experiments have been 
conducted in which the condensate was released into the 
waveguide using a method described previousl}i^. By 
carefully choosing initial conditions and s-wave scattering 
lengths, soliton trains were produced via a modulational 
instability mechanism in the single beam waveguide. Fig¬ 
ure [^sup) shows that the imaging system is clearly able 
to resolve large soliton trains. 

Finally, it is highly unlikely that the symmetric 
monopole mode of the in-trap oscillation extracted by 
PC A in Figure l(sup) would arise from a multicomponent 
soliton. The reasoning is twofold: firstly, it is unlikely 
that the breathing oscillation of each individual soliton 
in the strongly trapped direction would be in phase; sec¬ 
ondly, it is unlikely that the breathing oscillation of each 
individual soliton would have the same frequency as the 
symmetric oscillation of the soliton train in the weakly 
trapped axis. 

Following a previous method and the same dimension¬ 
less parameterJ^, the Gross-Pitaevski energy functional 
was used to calculate the energy surface with respect to 
both the axial and radial size of the condensate for the 
cubic-quintic GP: 

= ^ 2 I « I 2/3 

27 ^ 2 67 I 24 1357r^7f74 

(3) 

where /3 = N‘^'Re[gz]m/ In figure ^sup) the 

energy surface is plotted for = — 20 ao and the system 
parameters given in the main text. With Re[ 5 f 3 ] = 0, no 
stable point is present and the widths would tend towards 
collapse. However, for Re[^ 3 ] = 1.23 x 10“^^ficm^s“^ a 
stable point is clearly present. The inclusion of a positive 
Re[^ 3 ] stabilises the condensate at large negative scatter¬ 
ing lengthJ^, consistent with experimental observation 
of breathers in this regime. 
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FIG. 6 . (Sup) Energy surface plots of equation]^ sup) against 
radial and axial width at as — — 20 ao for both a) ^3 = 0 
and b) Re[^ 3 ] = 1.23 x (b). Parameters are: 

{N=1.5 X 10^, ujz = X 5.88 Hz, ujr = 27v x 77 Hz}. For 
R 6 [ 5 ' 3 ] = 0 the energy surface tends towards collapse, for non 
zero Re[^ 3 ] a stable point is shown. 
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FIG. 7. (Sup) In further experiments, soliton trains were pro¬ 
duced via modulational instability in a single beam optical 
waveguide. Total atom number in these images is 1.5 x 10^, 
and expansion time from the trap is 20ms. The imaging sys¬ 
tem is clearly able to resolve multiple solitons. 

experimental data in Figure |^. The phase shift between 
the axial and radial oscillations of 0.257r agrees with the 
measured shift of (0.24 ± 0.06) tt. The same expansion 
with Re[ 5 f 3]=0 does not agree with the experimental 
data (Figure]^). 


V. Phase shift arising from ballistic expansion 

As the cloud is initially released from the optical 
trap at high density, repulsive three body interactions 
during the dynamics of the expansion change the phase 
of the measured cloud widths. It is thus important to 
include when simulating ballistic expansion of the 
cold cloud prior to imaging. Figure |^sup) shows a 
simulation of the quadrupole oscillation for the same 
conditions as the experimental data in Figure 2 of 
the main text. The simulated expansion with = 0 
and Re[^ 3 ]= 1.23 x 10“^^ficm^s“^ is plotted with the 
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(a) In trap 



FIG. 8 . (Sup)Axial and radial widths of a simulated BEG 
with N = 1.5 X 10^, undergoing quadrupole oscillations at 
— 13ao inside a cigar shaped trap {ujz = 27rx5.88, ujr = 27rx77 
Hz), a) The in trap widths of the condensate showing the out 
of phase oscillations characteristic of the quadrupole mode. 

b) The widths of the condensate after undergoing 20ms of 
expansion with = 0 and Re[^ 3]=0 at each point in time. 

c) The expanded widths with Re[^ 3 ]= 1.23 x 10“^^/icm®s“^ 
plotted with the experimental data 
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